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15.  :  "rrum  2C0  vVcrasi 

Presented  here  in  this  Scientific  Report  Number  1  axe  the  preliminary  results  of  the  first 
phase  of  our  two-year  investigation  into  the  regional  wave  attenuation  in  Eurasia.  Our 
first  six-month  effort  ha.s  focused  upon  the  Lg  waves,  as  recorded  by  the  IRIS  and  CDSN 
stations.  A  time-domain  method  for  measuring  the  Lg  attenuation  has  been  developed 
and  successfully  tested  in  Eurasia,  yielding  stable  results  which  are  directly  compatible 
with  Nuttli’s  method  for  seismic  magnitude/yield  estimation.  While  continuing  to  refine 
our  Lg  results  through  improvements  in  path  coverage  and  methodological  upgrading,  our 
future  efforts  will  be  directed  at  gaining  a  better  understanding  of  the  P„  attenuation  in 
the  region. 
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SUMMARY 

The  pioneering  work  on  Lg  magnitude,  initiated  by  Nuttli  in  the  early  1970s,  has  been 
seen  blossoming  in  recent  years  into  notable  advancements  in  the  use  of  Lg  amplitudes  for 
precise  explosion  yield  estimation  (Nuttli,  1986,  1988),  propelling  intense  search  for  a 
reliable  LgQ  measurement  method.  Numerous  recent  publications  on  this  subject  involve 
mostly  spectral  analysis  of  Lg  wave  trains  of  varying,  somewhat  arbitrarily  chosen  group- 
velocity  windows  or  fixed- length  time  windows  {e.g..  Hasegawa,  1985;  Chun  et  al.,  1987; 
Baumgardt  and  Ziegler,  1988;  Sereno  et  a/.,  1988;  Sereno,  1990).  This  program  seeks  to 
determine,  eimong  others,  reliable  LgQ  in  time  domain  for  Eurasia,  rendering  the  result 
more  directly  compatible  with  Nuttli’s  method  of  explosion  yield  estimation. 

We  report  here  a  new  time-domain  method  for  measuring  Lg  attenuation.  Termed 
source  pair/receiver  pair  (SPRP)  method,  the  salient  features  of  this  approach  are;  a) 
effective  decontamination  of  recording  site  effects,  and  instrument  response  error;  and 
b)  effective  removal  of  source  excitation  function.  We  also  report  our  preliminary  Lg 
attenuation  measurement  results  for  Eurasia  based  upon  the  data  from  the  IRIS  and  CDSN 
stations.  Short  IRIS  and  CDSN  network  recording  histories,  high  event  trigger  threshold 
for  the  latter,  and  sparse  network  station  coverage  are  factors  which  limit  the  spatial 
resolution  of  our  regional  Lg  attenuation  measurement.  Nevertheless,  the  application  of 
our  new  methodology  has  yielded  surprisingly  robust  regional  results.  The  stability  of  the 
results  implies  higher  attainable  spatial  resolution  (fine-scale  regionalization)  as  the  Lg 
path  coverage  improves.  The  path  coverage  improvement  will  come  about  as  a  result  of 
our  continuing  acquisition  of  new  IRIS  and  CDSN  data  during  the  remaining  18-months 
period  of  our  two-year  research  project,  zis  well  as  the  use  of  additional  Chinese  seismic 
data. 
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1.  INTRODUCTION 

In  the  axea  of  regional  verification  of  nuclear  test  ban  treaty  compliance  the  DARPA 
research  program  has  focussed  sharply  in  recent  years  upon  the  use  of  Lg  phase  for  both 
source  discrimination  and  explosion  yield  estimation.  We  note  that  LgQ  estimates  derived 
from  the  observed  spectral  amplitudes  abound  in  the  published  seismological  literature. 
These  estimates  axe  most  useful  for  the  spectral-ratio  types  of  explosion/earthquahe  source 
discrimination  {e.g.,  Lg/Pn).  It  appears,  however,  that  LgQ  estimate  obtained  in  the 
time  domain,  talcing  into  consideration  of  the  Airy  phase  nature  of  the  Lg  wave  {e.g., 
Nuttli,  1973;  Campillo  et  ai,  1984),  provides  a  more  pertinent  attenuation  correction 
term  in  the  magnitude/yield  calculations  (Francis  Wu,  personal  communication,  1992). 
This  is  because  the  Nuttli’s  magnitude/yield  calculations  are  based  upon  time-domain  Lg 
amplitude  measurements  made  at  around  1-sec  period  rather  than  spectral  ampUtudes 
derived  over  finite  time  windows  (e.g.,  17-sec  and  34-sec  window  lengths  were  used  in 
Hasegawa,  1985;  3. 6-3.0  km/sec  group  velocity  window  was  used  in  Sereno,  1990). 

Here  we  describe  a  method  for  measuring  Lg  attenuation.  We  demonstrate  the  relia¬ 
bility  as  well  as  the  flexibility  of  the  method  by  applying  it  to  73  Lg  recordings  produced 
by  10  large-magnitude  earthquakes  and  registered  at  9  IRIS  and  CDSN  stations  (Table  1 
and  Figure  1).  The  LgQ  results  presented  in  this  report  are  of  preliminary  nature.  More 
refined  results  on  regional  wave  attenuation,  including  both  Lg  and  P„  Q,  will  be  presented 
in  future  reports. 


2.  A  TIME-DOMAIN  METHOD  FOR  MEASURING  Lg  ATTENUATION 

In  what  follows  we  describe  the  source  pair /receiver  pair  (SPRP)  method  for  measur¬ 
ing  regional  Lg  wave  attenuation. 

We  denote  by  Aij  the  time-domain  Lg  amplitude  due  to  the  source  and  recorded  at 
the  station.  We  assume,  following  Nuttli  (1973)  and  inserting  instrument  amplification 
Ij  with  its  attendant  error,  site  amplification  Sj,  and  source  radiation  pattern  Ri,  that 
Aij  can  be  parameterized  as 


Ai,,  =  (1) 

where  Ai  is  the  source  excitation,  dij  the  epicentral  distance  in  ki^jmeters  between 
the  source  and  station,  A  is  the  same  epicentral  distance  but  measured  in  degrees, 
and  7  the  coefficient  of  anelastic  attenuation.  This  coefficient  is  related  to  the  quality 
factor  Qif)  by  the  equation 
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7(/)  = 


UQifY 


(2) 


where  U  is  the  Lg  group  velocity. 

Consider  the  diagram  in  Figure  2  where  sources  1  and  2,  are  aligned  with  stations  1 
and  2  along  a  great-circle  path  (shown  as  a  solid,  horizontal  line).  This  configuration  is 
reminiscent  of  the  reversed  two-station  method  (RTSM)  deployed  by  Chun  et  al.  (1987), 
and  Zhu  et  al.  (1991)  in  their  regional  phase  attenuation  studies  in  eastern  Canada  -  except 
for  the  fact  that  the  two  stations  now  straddle  the  source  pair  rather  than  the  other  way 
around.  It  can  be  readily  shown  that 


/  d.2.1  ■^1,2  \  / <^2,1  di,2 
V-^1,1  ^2,2/  Vdi.i  ^2,2 


1/3 


/ sinA2,i  sinAi,2 
\5inAi,i  sinA2,2 


1/2 


=  e — 


(3) 


where  D\  —  <^2,1  —  and  D2  —  2  —  (^2,2- 

Note  that  in  talcing  the  Lg  amplitude  ratio  A2,i/Ai^i,  the  site  amplification  5i  and 
instrument  response  error  (if  any)  at  station  1  both  vanish  through  self-cancellations. 
Likewise,  all  source  excitation  and  radiation  pattern  terms  are  cancelled  when  the  product 
is  taken  of  the  amplitude  ratios  A2,ifAi^\  and  Ai,2/-A2,2-  As  indicated  schematically  in 
Figure  2,  a  partial  Lg  blockage  occurring  near  a  station  (KIV  in  Figure  1  for  instance) 
does  not  significantly  affect  our  Lg  attenuation  meeisurements. 

Calling  the  lefthand  side  of  equation  (3)  Y  and  taking  the  logarithm  on  both  sides 
lead  to 


logY  =  —^logeD, 


(4) 


where  D  =  Di  +  D2- 

For  simplicity,  we  now  make  the  usucil  assumptions  that  the  source  radiation  pattern 
and  the  site  amplification  are  both  azimuth  independent  for  the  Lg  wave  (e.^.,  Sereno, 
1990;  Bennett  et  al.,  1990).  Such  Eissumptions  axe  reasonable  since  the  Lg,  being  the 
superpositions  of  many  different  higher- mode  surface  waves,  is  presumably  more  capable 
of  “smoothing”  out  the  azimuth-dependent  effects  at  the  source  eind  receiver  ends  than  do 
the  body  waves.  Under  these  assumptions,  the  sources  and  receivers  axe  no  longer  required 
to  be  co-lineax  (see  dashed  propagation  paths  in  Figure  2).  Instead  of  a  variation  of  the 
RTSM  method,  we  call  this  generalized  version  source  pmr/receiver  pair  (SPRP)  method. 
Equation  (4),  which  describes  a  straight  line  with  zero  intercept  on  the  Y-axis,  remains 
applicable. 
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3.  PRELIMINARY  Lg  ATTENUATION  RESULTS  IN  EURASIA 


Figure  3  shows  four  seismograms  illustrative  of  the  quality  and  the  characteristics  of 
the  Lg  signatures  we  have  been  analyzing.  Two  are  IRIS  (GAR  and  OBN)  recordings 
and  rest  CDSN  (HIA  and  BJI)  recordings.  All  are  vertical-component,  broadband  data 
sampled  at  20  points/sec.  Only  distant  (A  from  2543  to  4964  km)  recordings  are  sampled 
here  to  show  that  the  Lg  waves  do  propagate  to  fax-regional  distances. 

The  actual  Lg  amplitudes  which  we  use  in  equations  (3)  and  (4),  are  time- 
domain  measurements  made  on  narrowband-filtered  (0.375  -  0.625  Hz)  records  (Figure  3). 
The  center  frequency  for  the  filter  is  0.5  Hz.  The  choice  of  this  particular  center  frequency 
is  based  upon  our  observation  that  the  signal-to-noise  ratio  (SNR)  at  0.5  Hz  is  good  for  the 
selected  propagation  paths,  long  and  short,  while  the  SNR  at  1-Hz  frequency  and  higher 
is  only  good  for  the  relatively  short  paths.  The  represent  the  maximum  zero-to-peak 
(trough)  Lg  amplitudes  taken  in  the  group  velocity  window  3.5  to  3.2  km  per  second  (see 
arrow  pairs  in  Figure  3). 

Shown  in  Figure  4a  is  a  plot  of  measured  logY  values  versus  D.  All  73  paths  are  used 
in  this  plot.  The  slope  of  the  fitted  least-squares  line,  as  discussed  previously,  is  -loge^. 
Because  of  our  concerns  for  the  potential  presence  of  systematic  bias  in  the  observed  logY, 
we  choose  not  to  fix  the  Y-intercept  at  the  origin,  as  implied  by  equation  (4).  In  other 
words,  we  effectively  allow  the  Y-intercept  to  be  a  free  parameter  of  the  least-squares  fit. 

The  slope  of  the  least-squares  line  shown  in  the  figure  gives  a  \’alue  of 
(1.129±0.062)xl0“^  for  7,  which  translates  to  a  Q  value  of  398,  2issuming  an  Lg  group 
velocity  of  3.5  km/sec.  Note  that  our  derived  Y-intercept  is  close  to  zero,  as  it  should. 

The  analysis  which  led  to  the  results  shown  in  Figure  4a  is  then  repeated  with  eight 
fewer  propagation  paths,  each  having  crossed  a  known,  near-station  partial  Lg  blockage. 
For  the  Lg  paths  crossing  the  northern  and  southern  Ural  Mountain  (Figure  1),  the  effects 
of  Lg  blockage  are  not  evident,  as  seen  on  our  filtered  seismograms.  The  deleted  ones  are: 

a)  five  paths  which  traverse  the  Caspian  Sea  to  reach  the  KIV  station  (Figure.  1);  and 

b)  three  Lg  paths  which  traverse  northern  Tibet  to  reach  the  KMI  station.  The  purpose 
here  is  to  examine  the  effects  of  near-station,  partial  Lg  blockages  on  the  measured  Lg 
attenuation. 

Figure  4b  shows  the  plot  of  logY  versus  D  with  the  eight  potentially  problematic 
paths  removed  from  Figure  1.  The  slope  of  the  least-squares  line  now  gives  a  value  of 
(1.045±0.078)xl0“^  for  7,  which  translates  to  a  Q  value  of  429.  The  result  given  in 
Figure  4b  is  our  preferred  Lg  attenuation  estimate  at  0.5  Hz. 

The  wide  epicentral  distance  span  (from  about  500  to  5.000  km),  necessitated  by  the 
large  inter-station  distances  and  non-ideal  epicentral  locations  of  the  earthquakes,  renders 
it  impossible  to  study  the  frequency  dependence  of  the  LgQ  iQ{f)  =  Qo/“).  since  the 
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high-frequency  signal  energies  axe  attenuated  along  the  paths  more  rapidly  than  their  low- 
frequency  counterparts.  Upon  reaching  a  fcir-regional  station,  only  the  low-frequency  part 
of  the  Lg  signal  energy  survives.  To  estimate  Qo  (Q  at  1  Hz),  it  is  then  necessary  to  know 
a.  Our  literature  search  reveals  that  in  tectonic  regions  Qo  is  small  and  a  is  large,  but  that 
the  opposite  is  true  in  stable  regions.  For  example,  the  Lg  Qo  in  the  Great  Basin  is  206 
and  the  a  0.68  (Chavez  and  Priestley,  1987);  the  corresponding  values  are  1,100  and  0.19 
in  eastern  Canada  (Chun  et  al.,  1987).  The  a  estimate  by  Given  et  al.  (1990)  for  eastern 
Kazakhstan  is  0.4.  According  to  Mitchell  (private  communication,  1992),  the  a  value  for 
Eurasia,  just  north  of  Tibet,  could  be  as  high  as  0.6.  Assuming  an  intermediate  o  value  of 
0.5  for  the  study  region,  our  preferred  result  (LgQ{0.5Hz)  seen  in  Figure  4b  would  then 
imply  a  Qo  value  of  607,  which  is  intermediate  between  206  for  the  Great  Basin  (Chavez 
and  Priestley,  1987)  and  1,100  for  the  Canadian  Shield  (Chun  et  al.,  1987). 


4.  CONCLUSIONS 


We  have  introduced  a  new  method  for  measuring  regional  Lg  attenuation.  Termed 
SPRP,  this  method  requires  no  great-circle  alignment  for  the  seismic  source  and  receiver 
pairs,  and  is  hence  fzu"  more  flexible  than  the  reversed  two-station  method  (RTSM)  of 
Chun  et  al.  (1087).  This  flexibility  allows  Lg  attenuation  measurements  to  be  made  using 
passive,  sparse  recordings  from  a  seismic  network  with  only  a  brief  recording  history. 

Like  the  RTSM,  the  SPRP  method  enables  effective  cancellations  of  the  source  ex¬ 
citation  function,  station  site  effects  and  is  capable  of  handling  near-station  partial  Lg 
blockage  to  some  extent.  We  have  successfully  applied  the  SPRP  method  to  broadband 
IRIS  and  CDSN  data.  As  shown  in  Figure  4a  and  4b,  the  source-  and  receiver-corrected 
observations  are  well-behaved  and  straightforward  to  interpret.  The  data  scatter  is  small 
by  normal  standards.  The  LgQ  at  O.o  Hz  is  429,  implying  a  value  of  607  if  a  is  taken 
to  be  0.5,  a  value  which  is  at  midpoint  between  Given’s  (1990)  value  of  0.4  and  Mitchell’s 
(personal  communication.  1992)  suggested  value  of  0.6. 

Compared  with  the  RTSM,  two  notable  drawbacks  of  the  SPRP  method  are:  a)  it 
does  not  permit  the  annihilation  of  the  source  radiation  pattern  effects;  and  b)  it  is  capa¬ 
ble  of  measuring  only  regional  attenuation,  not  path  attenuation.  We  wish  to  point  out, 
however,  that  for  Lg  waves,  which  consist  of  superpositions  of  higher-mode  surface  waves, 
the  radiation  pattern  effects  are  probably  not  very  important.  Also,  ais  more  data  accumu¬ 
late  with  time,  the  path  coverage  will  improve,  making  possible  fine-scale  Lg  attenuation 
mapping  and  consequently  the  derivation  of  effective  path  attenuation. 
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TABLE  1 


Event 

No. 

Date 

Origin  Time 
(Hr;Min:Sec) 

Latitude 

(N) 

Longitude 

(E) 

Depth 

(km) 

Magnitude 

(mb) 

1 

14  Jun  1990 

12:47:29 

47.87 

85.07 

58 

6.1 

2 

12  Nov  1990 

12:28:52 

42.96 

78.07 

19 

5.9 

3 

25  Feb  1991 

14:30:28 

40.39 

78.96 

21 

5.5 

4 

15  May  1990 

22:29:59 

36.11 

100.12 

14 

5.5 

5 

26  Apr  1990 

09:37:45 

36.24 

100.25 

10 

6.3 

6 

14  Jan  1990 

03:03:19 

37.82 

91.97 

12 

6.1 

7 

22  Sep  1989 

02:25:51 

31.58 

102.43 

15 

6.1 

8 

17  May  1989 

05:04:36 

57.09 

122.02 

31 

5.6 

9 

20  Apr  1989 

22:59:54 

57.17 

121.98 

26 

6.1 

10 

05  Mar  1990 

20:47:01 

36.91 

73.02 

12 

5.8 

Figure  1.  Lg  great-circle  paths  linking  the  epicentral  locations  of  the  ten  selected  earth- 
quaikes  (Table  1)  and  the  IRIS/CDSN  stations. 


PB  :  Partial  Blockage 

Figure  2.  Schematic  drawing  illustrating  the  source  pair/receiver  pair  (SPRP)  method 
for  measuring  Lg  attenuation.  Solid  horizontal  line  shows  an  ideal  great-circle  align¬ 
ment;  dashed  lines  show  the  “generalized”  configuration  which  is  actually  used  in  the 
study.  The  ideal  configuration  allows  czincellation  of  the  source  radiation  effects;  the 
generjilized  configuration  does  not. 
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Figure  4.  logF  (lefthand  side  of  equation  4)  plotted  as  a  function  of  D  (see  text).  The 
slope  of  the  least-squares  line  is  -logexj.  The  result  derived  from  all  73  paths  is  shown 
in  4a;  the  result  derived  from  65  selected  paths  (see  text)  is  shown  in  4b.  The  latter 
is  our  preferred  result. 
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